Dependence of spin torque switching probability on electric current 
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Dependence of the thermally assisted spin torque switching probability on the sweep electric 
current was investigated theoretically. The analytical expressions of the switching times for 6=1 
and b = 2 are derived based on the rate equation, where b is the exponent of the current term in 
the switching rate. The switching current is approximately proportional to the temperature T and 
the logarithm of the sweep rate v for both 6=1 and 6 = 2 in the experimentally performed ranges 
of T and v. Experiments in very low temperature range are required to determine the exponent 6. 



I. INTRODUCTION 

Spin random access memory (Spin RAM), which em- 
ploys spin torque switching P, 0] as the writing method, 
is one of the key spin-electronics applications for future 
nanotechnology. In Spin RAM application, a high ther- 
mal stability Ao of more than 40 is required to guarantee 
retention time longer than 10 years. Recently, Hayakawa 
et al. and Yakata et al. studied the thermal stabilities of 
antiferromagnetically and ferromagnetically coupled syn- 
thetic free layers respectively, and found higher thermal 
stabilities compared with that of a single free layer . 

The thermal stability is determined by the thermally 
assisted magnetization switching probability under the 
effect of the applied magnetic field or spin torque. The 
magnetic field or spin torque is assumed to be small 
compared to the anisotropy field Hk or critical current 
7 C at T = 0. In 1963 Brown derived the analytical 
expression of the switching probability which includes 
the effect of constant field -ff a ppi [6j which is given by 
P = l-exp{-/ iexp[-A (l-.f7appi/#K) 2 ]}, where f is 
the attempt frequency. The switching probability under 
the effect of constant current I (spin torque) was inde- 
pendently derived by Koch et al. and Li and Zhang [1] 
in 2004, and is given by P = 1 — exp{— fob exp[— Aq(1 — 
H app \/H K ) 2 (l - I/Ic)}}- It should be noted that the 
exponent of the current term, (1 — I/I c ), in this for- 
mula is unity. On the other hand, Suzuki et al. showed 
that the switching probability is given by P = 1 — 
exp{-/ texp[-A (l - H, ppl /H K ~ I/I c ) 2 }} [9], where 
the exponent of the current term is 2. Recently, we de- 
rived the switching probability formula which is given by 
P = 1 - exp{-/ i exp[-A (l - H appl /H K ) 2 (1 - I/I c ) 2 }} 
[1C| . which slightly different compared to the result of 
Suzuki et al. but agrees well with the results of Butler et 
al. [TTj] . The determination of the value of the exponent 6 
of the current term in the switching probability formula, 
P = 1 - exp{-/ t exp[-A (l - H appl /H K ) 2 (1 - I/I c ) b }}, 
remains unsettled, although it is very important for the 
determination of the thermal stability [Hj . 
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Recently we [T(| pointed out the reason why the values 
of the exponent 6 for the spin torque switching are dif- 
ferent between Ref. [7[ and Ref. [10( . In the case of the 
completely in-plane switching of the in-plane magnetized 
system, or in the case of the perpendicularly magnetized 
system, the switching can be described by only the angle 
from the easy axis, and the effect of spin torque can be re- 
garded as the enhancement of the applied field, as shown 
in Ref. [9T-fli"j|. Thus, the exponent b is reduced to 2, as 
is in the case for the field switching Q. However, the 
results of Refs. [tIB (b = 1) are supported by the exper- 
imental result |12j |. in which the mean switching current 
is approximately proportional to the temperature T and 
the logarithm of the current sweep rate log v. 

In this paper, we derived the analytical expressions of 
the switching times t sw for 6 = 1 and b — 2 due to the 
time dependent current I(t), and calculated the depen- 
dence of the switching current I(t sw ) on the temperature 
T and the sweep rate v. We found that I(t aw ) is approx- 
imately proportional to T and \ogv for both 6=1 and 
6 = 2 in the experimental ranges, and thus, the value 
of the exponent cannot be determined by the results of 
Ref. [l2j]- We also found that a low temperature exper- 
iment will clarify the value of the exponent 6. I(t sw ) is 
proportional to T if 6 = 1 while it is nonlinear to T if 
6 = 2. 

This paper is organized as follows. The analytical ex- 
pressions of the switching probability P and the switch- 
ing time t sw for 6=1 and 6 = 2 are derived in Sec. |H] 
and in Sec. IIII1 respectively. The dependences of the 
switching current I(t sw ) on the temperature T and the 
sweep rate v are discussed in Sec. IIVI Section [V] is the 
summary of this paper. 



II. SWITCHING PROBABILITY 

In this section, we derive the analytical expression of 
the magnetization switching probability due to the sweep 
current (6 = 1 and 6 = 2). We consider spin torque 
switching of the free layer magnetization in a ferro- 
magnetic (fixed layer) /nonmagnetic/ferromagnetic (free 
layer) trilayer. Both ferromagnetic layers are assumed to 
have the perpendicular anisotropy along the z axis. The 
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magnetization M of the free layer is assumed to point to 
the positive z direction at t — 0. Throughout this paper, 
the applied field H app \ is assumed to be zero for simplic- 
ity. From t = 0, the current / is applied to the ferro- 
magnetic film perpendicular to plane and induces torque 
on the magnetization M. The magnitude of the current 
I is assumed to be smaller than the critical current I c . 
With the help of the thermal activation, the magnetiza- 
tion can change its direction to the negative z direction. 
The switching probability of the magnetization, P, obeys 
the following equation [6]: 



1.0 r 



dt 



r(t) [1-P(t)] 



(1) 



where r(t) is the switching rate from the initial state 
(M || e z ) to the final state (M || — e 2 ). The explicit form 
of r(t) is given by 



r (t) = fo exp 



-A 1 



I{t) 



(2) 



Here b is the exponent of the current term (1 — I/I c )- 
Although the attempt frequency depends on I{t) in gen- 
eral, we approximate the value of fo as the attempt fre- 
quency at zero current. This approximation is applica- 
ble because the dependence of the switching rate r(t) on 
the current I(t) is mainly determined by the exponential 
term exp[— Ao(l — I/I c ) ]. The attempt frequency with 
zero current is given by fo — cry/fx \/ Ao/7r/ (1 + a 2 ), 
where a and 7 are the Gilbert damping constant and the 
gyromagnetic ratio, respectively [6j. 

Let us consider the thermally assisted magnetization 
switching due to time dependent current I(t). As is in 
experiments [H, Hi3 |. we assume that the strength of I(t) 
increases linearly with time, i.e., I(t) — vt, where v is the 
sweep rate. By integrating Eq. (fTJ) and using the initial 
condition P(0) — 0, P(t) is given by 



P(t) = 1-exp. 



/q/ c 
bvA 



i/b 




(3) 

where j(/3, z) = f^dtt l3 ~ 1 e~ t is the incomplete T func- 
tion. For b = 1 1 and b = 2 @, UuM, p 0) are , 
respectively, reduced to 



P^t) = 1-exp 



P 2 (t) = 1-exp 



vA 



exp 



-A 1 



vt 



(4) 



exp(-Ao)]}, 
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FIG. 1: The time evolutions of P(t) for b ■■ 
and b = 2 (P 2 in Eq. ©). 



1 (Pi in Eq. ©) 



where Erf(z) = (2/1/7F) J* die - ' 2 is the error function. 
Equations (j4} and ([5]) are the main results in this section. 

Figure [T] shows the time evolutions of Pi and P2 un- 
der the effect of the sweep current. The values of the 
parameters are taken to be a = 0.007, 7 = 1.764 x 10 7 
Hz/Oe, H K = 200 Oe, M = 995 emu/c.c, T = 300 
K, d = 2.0 nm and S — n x 80 x 35 nm 2 , where d 
and S are the thickness and cross section area of the 
ferromagnetic layer, respectively, which are typical val- 
ues for the Spin RAM cell consists of the transition 
ferromagnetic metals such as CoFeB 0, H|. The ther- 
mal stability is defined by Ao = M HkV / '(2/cbP) , where 
V = Sd is the volume. The critical current J c is given 
by I c = [2aeMSd/{hg)]{H K + 2ttM) [l|, where spin 
polarization g is taken to be 0.5. By using the above pa- 
rameters, the values of fo, A , and I c /S are given by 0.09 
GHz, 42, and 5.5 x 10 6 A/cm 2 , respectively. The sweep 
rate is taken to be v/S = 5.0 x 10 6 A/cm 2 s (v = 0.44 
mA/s), which is similar to the experimental values [12j . 

The switching probability P suddenly changes its value 
from to 1 at a certain time t sw , as shown in Fig. [TJ in 
which t sw satisfies I(t BW ) < I c . In next section, we derive 
the analytical expression of the switching time i sw . 



III. SWITCHING TIME 

In this section, we derive the analytical expression 
of the switching time. At t = i sw , dP/dt takes its 
maximum, as shown in Fig. [1] Thus, P(t) satisfies 
d 2 P/d£ 2 = at t — t sw . This condition can be rewritten 
as 



dr 
dt 



(6) 



By using the explicit form of r(t) (Eq. ©), we find that 



bvA 

T, 



i{t sw ) 



b-l 



/oexp 



-A 1 



I{t sw ) 



(7) 

In the case of b = 1, Eq. can be easily solved, and we 
find that 



I(ts 



1 : — lOET 



A vA 



(8) 
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FIG. 2: (a) The dependences of the switching times for 6=1 
and & = 2 (Eqs. © and (V2\i ) on the thermal stability A . 
(b) The relation between the thermal stabilities estimated by 
Eqs. (f9)> and (fT2)). 

By using / = vt, the switching time for 6 = 1 is given by 



1 , f Ic 

Ao 10g ^Au 



(9) 



which is identical to that obtained by Li and Zhang [8J 
for P = 1 - e" 1 (see Eq. (22) in Ref. 0). 

To obtain the explicit form of t S w for general 6, with the 
help of Eq. ((5J, let us consider the solution of I(t sw )/I c 
whose form is given by 

1/6 



1 



1 

An 



\ogX 



(10) 



where X is determined by Eq. ([7]). By substitut- 
ing Eq. (fl"0|) into Eq. (UJ), we find that X satisfies 

(logX) 1 - 1 ^ = 1/(CX), where C = bv^ /b / '(/ I C ). In 
terms of Y = logX, this equation can be expressed as 
Ye Y Y~ 1 / b — 1/C, and its solution is given by Y = 
[{b - l)/b]p\og[bC~ b ^ b '^/(b - 1)], where plog(z) is the 
product log (or Lambert W function) which satisfies 
plog(z) exp[plog(z)] = z. I(t sw )/I c is thus given by 



/(t sw ) 



= 1-. 



6-1 
~6Ao~ 



plog 



Zoic 



6/(6-1) 



6vA 



1/6 



1/6 



In the case of b = 2, the switching time t s 



is given by 



k 

v 







h 


2A plOg [ 



/do 



(12) 



Equations ([9]) and (fT2"|) are the main results in this sec- 
tion: these are the theoretical expressions of the switch- 
ing times t sw , or equivalently, the switching current 
I (t S w) = vt sw - They are valid for the sweep rate v which 
satisfies i sw > 0. By using the parameter values listed in 
the previous section, i S wi and i SW 2 are estimated to be 
0.71 s and 0.45 s, respectively, which have good agree- 
ment with Fig. [T] 

For a large z, plog(z) can be approximated to log(z). 
In this limit, Eq. (jlll) is reduced to 



I{t s 



1 - 



aV os 



6-1 



1-1/6 



M 

6wA 



1/6 



Neglecting a small correlation term [6/(6 — l)] 1_1//b , 
Eq. (fT3"l) is identical to the result of Garg [l4| 
on the mean switching current (/) / I c — 1 — 
{(l/A )log[/o/ c /(6wAo /b )]} 1 / b which was used in the 
analysis of the experimental result of Ref. [l2[ • 

Experimentalists [1, [l^] repeat the measurement of the 
switching due to sweep current or magnetic field and 
measure the switching time of each trial. By fitting the 
average of the switching time with Eqs. © or (jT2"j) . one 
can estimate the thermal stability Ao- Figure[2](a) shows 
the dependence of the switching times, Eqs. (j9]) and (|12[) . 
on the thermal stability Aq, in which Ao is regarded as 
the independent variable from J c , as is in experiments 
For example, when the switching time is about 0.8 
s, the thermal stability estimated by Eq. © is about 50 
while that estimated by Eq. (IT21 is about 200. Figure 
[2] (b) shows the relation between the thermal stabilities 
estimated by Eqs. © (Ai) and (fH?)) (A 2 ) which can be 
obtained by putting t sw i(Ao = Ai) = i sw2 (A = A 2 ) 
and is given by 



plog 



/o/c 



\ 



plog 



hie 




(14) 

As shown in Fig. [2] (b), the difference of the exponent 
6 leads to a significant underestimation of the thermal 
stability. 



IV. DEPENDENCE OF SWITCHING TIME ON 
TEMPERATURE AND SWEEP RATE 



Myers et al. experimentally studied the dependences of 
the switching current, I(t sw ), on the temperature T and 
the sweep rate v (see Fig. 1 (d) in Ref. [Li), and showed 
that /(t sw ) is approximately proportional to both T and 
log v. Reference [8[ argued that these results support the 
theory with 6 = 1 [7, 8j because such dependences are 
expected from Eq. (|8]). In this section, we study the 
dependences of I(t S w) for 6 = 1 and 6 = 2 on T and v in 
and out of the experimental ranges. 

The dependence of I(t sw ) on the temperature can be 
taken into account through the thermal stability Ao = 
M HkV I '(2k-QT) . For simplicity, we neglect the depen- 
dence of the attempt frequency /o oc a/Ao on the temper- 
ature T, as done in Ref. [8]. The solid (red) and dotted 
!/ b (blue) lines in Fig. [3] show the dependences of I(t S w)/I c 

, on T for 6=1 and 6 = 2, respectively. The temperature 

range (100-300 K) is sufficient to consider the I(t sw ) in 
(13) the experimental range (180-220 K in Ref. [3). Even in 
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FIG. 3: The solid (red) and dotted (blue) lines are the de- 
pendences of I(t s -w)/Ic on the temperature T for b — 1 and 
b — 2, respectively. The dashed (orange) and dashed-dotted 
lines are the dependences of I(t avv )/I c on the sweep rate v for 
6=1 and 6 = 2, respectively. 

(b) 
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FIG. 4: The dependences of I(t av ,)/I c on (a) temperature 
T and (b) the sweep rate v out of the experimental range of 
Ref. Qj]. 

such wide range of the temperature, both I(t sw )/I c for 
6 = 1 and 6 = 2 are approximately proportional to the 
temperature T. The dependences of I(t sw )/I c on u for 
6=1 and 6 = 2 are also plotted by the dashed (orange) 
and dashed-dotted (green) lines in Fig. [3J respectively. 
The temperature is fixed to T — 300 K. The range of 
v/S, from 0.05 to 5.0 xlO 6 A/cm 2 s (from 0.0044 to 0.44 
mA/s), is very similar to that studied in experiments 
(from 0.01 to 1.0 mA/s). As shown in Fig. [31 both 
I(t S w)/Ic for 6=1 and 6 = 2 are approximately propor- 
tional to logu. Summarizing the above results, the value 



of the exponent 6 (6 = 1 or 6 = 2) cannot be determined 
by the dependences of I(t aw )/I c on the temperature T 
and the sweep rate v in the experimental range. 

We also study the dependences of I(t sw )/Lon T and 
log v out of the experimental ranges of Ref. [12| . Figure 
2] (a) shows the dependence of I(t sw )/I c on T in the low 
temperature region, < T < 50 K. In this range, we can 
see that I(t sw )/I c is approximately proportional to T for 
6=1 while I(t S vi)/I c for 6 = 2 is not. Thus, to deter- 
mine the exponent 6, experiments should be performed 
in such very low temperature region. The dependences 
of J(t sw )/J c on logv for both 6=1 and 6 = 2 are ap- 
proximately proportional to log v even in the large sweep 
rate region up to v/S — 500 x 10 6 A/cm 2 s, as shown in 
Fig. [4] (b), and thus, it is difficult to determine 6 from 
these dependences. 

V. SUMMARY 

In summary, we studied the dependence of the spin 
torque switching probability on the sweep current. The 
analytical expressions of the switching time for both 
6=1 and 6 = 2 are derived, where 6 is the exponent of 
the current term in the switching rate. We showed that 
I(t s -w)/Ic for 6=1 and 6 = 2 are approximately pro- 
portional to the temperature T and the logarithm of the 
sweep rate v in the experimental ranges of T and v, and 
thus, the value of the exponent 6 cannot be determined 
by the dependences of I(t sw )/I c on T and log?;. We also 
showed that a low temperature experiment is required to 
determine the exponent 6. I(t svr ) is proportional to T if 
6=1 while it is nonlinear to T if 6 = 2. 
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